WEIGHTED ENUMERATION OF SPANNING 
SUBGRAPHS WITH DEGREE CONSTRAINTS 



DAVID G. WAGNER 

Abstract. The Heilmaim-Lieb Theorem on (univariate) match- 
ing polynomials states that the polynomial ^2 k irik(G)y k has only 
real nonpositive zeros, in which mj.(G) is the number of fc-edge 
matchings of a graph G. There is a stronger multivariate version 
of this theorem. We provide a general method by which "theo- 
rems of Heilmann-Lieb type" can be proved for a wide variety of 
polynomials attached to the graph G. These polynomials are mul- 
tivariate generating functions for spanning subgraphs of G with 
certain weights and constraints imposed, and the theorems specify 
regions in which these polynomials are nonvanishing. Such the- 
orems have consequences for the absence of phase transitions in 
certain probabilistic models for spanning subgraphs of G. 



1. Introduction. 

Let G = (V, E) be a finite graph, possibly with loops or multiple 
edges. For each natural number k e N, let m^G) denote the number 
of A;-edge matchings in G. The univariate Heilmann-Lieb Theorem [I] 
states that all zeros of the polynomial fi(G;y) = J2k m k(G)y k lie on 
the negative real axis. A stronger multivariate version has variables 
x = {x v : u6^}, one for each vertex, and concerns the polynomial 

Ji(G;x) = J2* deg{M) 

M 

in which the sum is over all matchings M of G, deg(M) : V — > N is the 
degree function of M, and for any / : V — > N 



n 
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The multivariate Heilmann-Lieb Theorem [lj states that if | &rg(x v )\ < 
7r/2 for all v G V then Jl(G; x) 7^ 0. One sees that this implies the 
univariate version by means of the relation 

^(G;y) = jl(G;y 1 / 2 l) 

(which follows from the Handshake Lemma). 

The purpose of this paper is to apply some standard results from the 
analytic theory of complex polynomials to provide a general method by 
which "theorems of Heilmann-Lieb type" can easily be deduced. The 
multivariate Heilmann-Lieb Theorem itself appears as the simplest - 
and prototypical - special case of the method. Other direct applications 
provide multivariate extensions of previous results of the author [TO] , 
and of results of Ruelle [H [9] . A variety of new results also appear as 
natural special cases. 

In the remainder of this Introduction we describe the general combi- 
natorial situation we will consider. In Section 2 we gather the necessary 
results from the analytic theory of complex polynomials. In Section 3 
we state and prove the main theorem of the paper. Section 4 illustrates 
this result with several applications, including the previously known ex- 
amples mentioned above. In Section 5 we explain an interpretation of 
the polynomials we consider as partition functions, by analogy with 
the Boltzmann-Gibbs formalism in statistical mechanics. Results like 
those in Section 4 imply that when the thermodynamic limit of the 
free energy exists it must be analytic in certain regions of the complex 
plane. As noted by Lee and Yang [5J [TT] , this has implications for the 
absence of phase transitions in these models (which enumerate span- 
ning subgraphs subject to certain weights and constraints). A more 
thorough investigation of the phase structure of these models would be 
very interesting, but must be left for a later paper. 

It is a pleasure to thank my good friend Alan Sokal, the anonymous 
referee, whose detailed positive criticism of an earlier form of this pa- 
per prompted me to rewrite it completely. The result is, I think, much 
improved. 

The general framework we consider is that of a finite graph G = 
(V, E) (possibly with loops or multiple edges) and a set of weights 
A = {A e : e G E} on the edges of G. These weights can for some 
purposes be considered as indeterminates, but will usually be taken to 
be complex numbers, and often will be nonnegative real numbers. (In 
combinatorial applications it is most natural to set all the edge-weights 
equal to one.) The starting point for the theory is the elementary 
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identity 

(1.1) Q(G, A; x) = 11 (1 + \ e x v x w ) = ^x dc6(if) . 

veweE HCE 

In this formula, the product is over the set of all edges e G E, and 
the notation vew indicates that the ends of e are the vertices v and w 
(note that v — w is possible). The sum is over the set of all spanning 
subgraphs (V, H) of G, each of which is determined by its edge-set 
HCE. As above deg(H) : V — > N is the degree function of if, and 
we use the shorthand notations 

eeH 

and 

This fi((j, A; x) is a relatively structureless object, since it sums over 
all spanning subgraphs without preference. On the other hand, the 
product formula allows one to make very precise statements about its 
zero-set (as a subset of C v ). To make use of this, we introduce a 
sequence of activities at each vertex v G V: 

(1.2) u<"> = (i4 w \u< w \ ...,«?>) (d = deg(G») 

which can be any complex numbers (usually taken to be nonnegative 
reals). With these activities specified, a spanning subgraph HCE 
will be given the weight 

( L3 ) U dcg(H) = J] U dc g (H,v) 

vev 

and we will consider the correspondingly weighted version of Q(G, A; x): 
(1.4) Z(G, A, u;x)=£ A H u dc#) x d ^. 

HCE 

For example, if at every vertex we take u — u± — 1 and -u^ = for all 
fc > 2, then 

J 1 if if is a matching, 
u de g(i/) - | o otherwise, 

and A, u; x) is an edge- weighted version of the multivariate match- 
ing polynomial J1(G; x) above. 

The strategy in what follows is to begin with information about the 
zero-set of Q(G, A; x) and to impose conditions on the vertex activities 
u^) that are sufficient to imply similar information about the zero-set 
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of Z(G, X, u;x). To realize this plan, we need a few results from the 
analytic theory of complex polynomials. 

2. Complex polynomials. 

The technique we use is known as Schur-Szego composition. We do 
not make use of the most general possible result, but for thoroughness 
of exposition we derive what is needed from the Grace-Szego- Walsh 
Coincidence Theorem. For a more complete treatment see Sections 15 
and 16 of Marden [6] and Chapters 3 and 5 of Rahman and Schmeisser 

m- 

Let F(z) be a polynomial in complex variables z := {z v : v £ V}. 
For a subset A C C, we say that F is A-nonvanishing if either F = 0, 
or z v £ A for all v £ V implies that F(z) ^ 0. In the case that F ^ 
we say that F is strictly .A-nonvanishing. 

Lemma 2.1. Le£ .A 6e nonempty, connected and open. Let F n (z) be 
a sequence of strictly A-nonvanishing polynomials indexed by positive 
integers, and assume that the limit F(z) = linv^oo F n (z) exists. Then 
F is A-nonvanishing. 

Proof. Each F n is analytic and strictly nonvanishing on the subset A v 
of C y . Since these functions are polynomials, the convergence to F is 
uniform on compact subsets of C v . By Hurwitz's Theorem (Theorem 
1.3.8 of [7j), either F is identically zero or F is nonvanishing on A v as 
well. ' " □ 

Lemma 2.2. Let A be nonempty, connected and open. Let F(z) be 
an A-nonvanishing polynomial, and let w £ V . If z w is fixed at a 
complex value £ in the closure of A, then the resulting polynomial in 
the variables {z v : v £ V \ {w}} is A-nonvanishing. 

Proof. The result is trivial if F = 0, so assume instead that F is strictly 
.A-nonvanishing. Let (£„ : n = 1, 2, ...) be a sequence with each £ n £ .A 
such that lim^oo = £. Note that for all n > 1 the specialization 
= results in a polynomial F n that is strictly .A-nonvanishing in 
the variables {z v : v £ V \ {w}}. The sequence (F n : n > 1) satisfies 
the hypothesis of Lemma 2.1, from which the result follows. □ 

We are concerned mostly with the following open subsets of C. 
• For < 9 < it, the open sector 

(2.1) S[6] = {zeC: z ^ and | arg(z)| < 6} 

centered on the positive real axis. (For z ^ we use the value of the 
argument in the range — ir < arg(z) < ir.) 
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• For k > 0, the open interior of a disk 

(2.2) kT> := {z G C : |z| < re}. 

• Also for re > 0, the open exterior of a disk 

(2.3) re£ := {z G C : \z\ > k}. 

When re = 1 we more simply write just D and £. 

A circular region in C is a proper subset that is either open or closed 
and is bounded by either a circle or a straight line. A polynomial 
F(z) = F(z\, Zd) is multiaffine if each variable occurs at most to 
the first power. The polynomial F{z) is symmetric if it is invariant 
under every permutation of the variables. The elementary symmetric 
functions of the variables z = (zi, z&) are 

(2.4) ej-(z) = ^ ZfcZi,,---^. 

I<ii<i2<---<ij<d 

A multiaffine symmetric polynomial F(zi, ...,^) is thus a linear com- 
bination of the elementary symmetric functions ej(z) for < j < d. 

Proposition 2.3 (Grace-Szego- Walsh). Let F(z\, z^) be a multi- 
affine symmetric polynomial, and let A be a circular region. Assume 
that either A is convex or the degree of F is d. Then, for any values 
Ci> •••> (d £ A there exists a value ( G A such that 

F(Ci,...,Cd = nc,...,C). 

For a proof in the case that degF = d, see Theorem 15.4 of [B] or 
Theorem 3.4.1b of [7]. The theorem also holds when degF < d with 
the additional hypothesis that A is convex, as explained in Theorem 
2.12 of pQ. 

For an elaboration of the ideas of Proposition 2.4, see Lemma 5.5.4 
and Theorem 5.5.5 of [7J. 

Proposition 2.4 (Schur-Szego). Let P(z) = Y2j c j z ^ an d K( z ) = 
Sj=o { t j) u j z '' be polynomials in one complex variable z, with degP < d, 
and let Q{z) = UjCjZ 3 . 

(a) For any < a < tc/2, if P(z) is S[7r /2}-nonvanishing and K(z) is 
§[n — a]-nonvanishing, then Q(z) is §[ti/2 — a]-nonvanishing. 

(b) For any k > and p > 0, if P(x) is pT>-nonvanishing and K{z) is 
kD -nonvanishing , then Q(z) is npT>-nonvanishing. 

(c) For any k > and p > 0, if P(x) is pZ-nonvanishing and K{z) is 
nt-nonvanishing and degK = d, then Q(z) is upE-nonvanishing. 
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Proof. The conclusions are trivial if Q = 0, so we may assume that 

We begin by proving part (a) in the case that K(0) ^ 0. In this case 
we have 

d 

(2.5) K(z)=cY[(l + 6 iZ ) 

i=i 

for some complex numbers C ^ and 9±, d such that either 9,i = 
or | arg(#j)| < a for each 1 < % < d. Consider the d-th polarization of 
P(z): this is the multiaffine symmetric polynomial -P(z) = P(zi, Zd) 
obtained from P(z) by replacing each monomial z^ by the normalized 

j-th elementary symmetric function (^) e_,-(z). Since degP < d, it 
follows that 

(2.6) P(z,z,...,z)=P(z) 

as polynomials in z. Since P(z) is S[7r/2]-nonvanishing and §[vr/2] is a 
circular region, it follows from (2.6) and Proposition 2.4 that P(z) is 
also §[7r/2]-nonvanishing. Now, consider complex numbers (\,...,(d £ 
S[tt/2 - a]. For each 1 < i < d, either ^0 = for all Q e S[tt/2 - a] 
or |arg(0j£j)| < 7r/2 for all Q e S[7r/2 — a]. From Lemma 2.2, it 
follows that if P(0iZi, ^^d) ^ then P(0iCi, OdQ ^ for every 
choice of Ci, ---Xd £ S[7r/2 — a]. That is, it follows that P(9iz 1 , ...,6 d Zd) 
is S[7r/2 — a]-nonvanishing. A short calculation using the fact that 
(f)Uj = Cej(6i, ...,9 d ) verifies that 

(2.7) Q(z) = CP(0 1 z,...,6 d z), 

and therefore Q(z) is §>[tt/2 — a]-nonvanishing, as desired. 

To handle the case in which K(0) = 0, let r be the multiplicity of 
as a root of K(z) and write 

d—r 

(2.8) K{z) = Cz r \[(l + 6 t z). 
For a positive integer N let 

d—r 

(2.9) K N (z) =CN~ r (l + Nz) r Y[(l + 6 t z). 

i=i 

and let Qn(z) be the polynomial in the conclusion constructed from 
P(z) and Kn(z). By the case we have done already, each Qn(z) is 
S[7r/2 — a]-nonvanishing. Taking the limit as iV — > oo, Lemma 2.1 
implies that Q(z) itself is also §[n/2 — a]-nonvanishing. 
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The proof of part (b) is similar. Since K(z) is /tD-nonvanishing we 
have K(0) ^ 0, and so we can write K{z) as in equation (2.5) with 
all \9i\ < 1/k. Again we consider the d-th polarization P(z) of P(z). 
Since P(z) is p!D-nonvanishing and pD is a circular region, Proposition 
2.3 and equation (2.6) imply that P(z) is p2)-nonvanishing. It follows 
that P(9iZi, ...,9dZd) is /tpD-nonvanishing, and from equation (2.7) we 
conclude that Q(z) is /tpD-nonvanishing, as desired. 

The proof of part (c) repeats the same pattern once more. Begin with 
K(z) expressed as in equation (2.8) - since K(z) is /t£-nonvanishing, 
each \9i\ > We work with the polynomials K^{z) defined in 

equation (2.9) with N > 1/k. Since P(z) is p£-nonvanishing and p£ 
is a circular region and degP = d, Proposition 2.3 and equation (2.6) 
imply that P(z) is p£-nonvanishing. It follows that 

P(9lZi, 9 d - r Z d -r, NZd-r+1, -, Nz d ) 

is /tp£-nonvanishing, and from equation (2.7) we conclude that Qn(z) 
is /tp£-nonvanishing. Taking the limit as N — > oo (using Lemma 2.1) 
we conclude that Q(z) is Kp£-nonvanishing, as desired. □ 

The polynomial Q(z) in the conclusion of Proposition 2.4 is the 
Schur-Szego composition of P(z) and K(z). 

3. The main result. 

Consider a graph G = (V, E) with complex edge weights A. We 
begin with some easy information about the zero-set of the polynomial 
£l(G, A;x) defined in equation (1.1). 

Proposition 3.1. Let G = (V, E) be a graph with complex edge weights 
A. 

(a) If X e > for each e G E then Q(G, A; x) is S[tt /2]-nonvanishing. 

— 1/2 

(b) If\X e \ < A max for each e G E then Q(G, A; x) is A ma x D-nonvanishing. 

— 1/2 

(c) If\X e \ > A m i n for each e G E then Q(G, A; x) is X mi ^ Z-nonvanishing. 

Proof. In each case, each factor 1 + X e z v z w in the product form for 
Q(G, A; z) is seen to be nonvanishing in the appropriate region, from 
which the result follows. □ 

Now assume that we also have a sequence of activities at each 
vertex v G V, as in equation (1.2). The information about these activ- 
ities that we will use is recorded in the set of key polynomials 

(3.1) = E (fW^' 

3=0 
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in which d = deg(G,v). There is one key polynomial for each vertex 
v G V. 

Theorem 3.2. Let G = (V, E) be a graph, with complex edge weights 
X, and with vertex activities u encoded by the key polynomials K v {z) 

(v e V). 

(a) Fix < a < n/2. // A e > for each e G E and K v (z) is 
§[n — a]-nonvanishing for each v G V , then Z(G, A, u; x) is S[n/2 — a\- 
nonvanishing. 

(b) Fix k > and A max > 0. If |A e | < A max for each e G E and K v {z) 

1 /2 

is nD-nonvanishing for each v G V , then Z(G, A, u; x) is (/t/Amax)2)- 
nonvanishing. 

(c) Fix k > and A min > 0. If |A e | > A min /or eac/i e E E and K v (z) 
is nZ-nonvanishing and deg K v (z) = deg(G, v) for each v G V, then 
Z(G,X,u;x) is (k/ \lli Q )8.-nonvanishing. 

Proof. Identify the vertices V with the numbers V = {l,2,...,n} ar- 
bitrarily. Define a sequence of polynomials -Fo(x), Fi(x),..., F n (x) as 
follows. -Fo(x) = Q(G, A;x), and for all 1 < v < n, F v (x) is the 
Schur-Szego composition of F^_i(x) regarded as a polynomial in the 
variable x v (the other variables being absorbed into the coefficients) 
with K v {x v ). One sees by induction that for < r < n: 

HCE \v=l J 

so that F n (x) = Z(G, A, u; x). 

We give the details to finish the proof of part (a) - the arguments 
for parts (b) and (c) are completely analogous. We prove by induction 
on 1 < v < n that if (Q : 1 < j < n) are complex numbers such that 

• | arg((j)| < 7r/2 — a for all 1 < j < v , and 

• | arg((j)| < 7r/2 for all v < j < n, 
then 

(3.3) F t) _ 1 ((' 1 , ( v -i, X v , Cv+l, Cn) 

is S[7r/2]-nonvanishing. The basis of induction follows from Proposition 
3.1(a) and Lemma 2.2. The induction step follows from Proposition 
2.4(a) and Lemma 2.2. Finally, from the statement that whenever all 
Q G §[n/2 — a], then F„_i(d, Cn-i, x n ) is S[7r/2]-nonvanishing, we 
conclude by one more application of Proposition 2.4(a) that -F n (x) is 
S[7r/2 — a]-nonvanishing, as desired. □ 

The univariate specialization of Theorem 3.2 is an important con- 
squence. 
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Corollary 3.3. Adopt the notation of Theorem 3.2. 

(a) Under the hypotheses of Theorem 3.2(a), 
Z(G, A, u; t/ 1 / 2 ]-) is S[ir — 2a]-nonvanishing. 

(b) Under the hypotheses of Theorem 3.2(b), 
Z(G, A, u; y 1 / 2 !) is (k 2 / X maiX )T)-nonvanishing. 

(c) Under the hypotheses of Theorem 3.2(c), 
Z(G, A, u; y x l 2 V) is (k 2 / \ min )£,-nonvanishing. 

4. Applications. 

Throughout this section, consider a graph G = (V, E) with complex 
edge weights A and vertex activities u encoded by the key polynomials 
K v (z) (veV). 

Example 4.1 (Heilmann-Lieb [1]). Assume that all edge weights are 
nonnegative reals, and that at each vertex uq — u\ — 1 and Uk = 
for all k > 2. The key polynomial at a vertex of degree d in G is 
K v (z) = 1 + dz, which is S[7r]-nonvanishing. Theorem 3.2(a) (with 
a = 0) implies that Z(G, A, u;x) is S[7r/2]-nonvanishing - this is the 
multivariate Heilmann-Lieb theorem. Corollary 3.3(a) implies that 
Z(G, A, u; 2/ 1 / 2 ].) is §[7r]-nonvanishing - this is the univariate Heilmann- 
Lieb theorem. 

Example 4.2 (Wagner [ID])- Assume that all edge weights are non- 
negative reals, and that two functions /, g : V — > N are given such that 
f{ v ) < g(v) < f{v) + 1 for each v G V. Fix the vertex activities to be 

(41) u («) = < [ 1 if f(v)<k<g(v), 

k 1 otherwise. 

As in Example 4.1, each key K v (z) is S[7r]-nonvanishing. Theorem 
3.2(a) (with a = 0) implies that Z(G, A, u;x) is S[7r/2]-nonvanishing 
- this result is new. Corollary 3.3(a) implies that Z(G, A, u; y^ 2 l) is 
S[7r]-nonvanishing - when A = 1 this is Theorem 3.3 of [TO] . 

Example 4.3 (Ruelle [U [9]). Assume that all edge weights are non- 
negative reals, and that at each vertex uq = ui = 1, u\ = u, and 
Uk = for all k > 2. The key polynomial at a vertex of degree <i in G 
is K v (z) = 1 + duz + (^j-? 2 - F° r d >2 } the zeros of this polynomial are 
at 

z± = (u ± ^u 2 - 2 + 2/d) . 

When u = 1, all the keys -fCj(z) are S[37r/4]-nonvanishing, and 
Theorem 3.2(a) (with a = vr/4) implies that Z(G, A, u; x) is §[7r/4]- 
nonvanishing - this is new. Corollary 3.3(a) implies that Z(G, A, u; y 1 ^ 2 !) 
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is S[7r/2]-nonvanishing - when A = 1 this is a slight weakening of 
Proposition 1 of [8]. 

If G has maximum degree A and u > a/2 — 2/A, then all the keys 
K v (z) are S[7r]-nonvanishing, and Theorem 3.2(a) (with a = 0) implies 
that Z(G, A, u; x) is S[7r/2]-nonvanishing - this result is new. Corollary 
3.3(a) implies that Z(G, A, u; |/ 1//2 1) is S[7r]-nonvanishing - when A = 1 
this is Proposition 2 of [8j. 

Ruelle's method produces more detailed information than ours, but 
only for particular choices of the vertex activities. A systematic exten- 
sion of his method that handles all the cases we consider would be very 
interesting. 

Example 4.4. Assume that all edge weights are nonnegative reals, 
and that two functions /, g : V — > N are given such that f(v) < g(v) < 
f(v) + 2 for each v G V. Fix the vertex activities as in equation 
(4.1). Then each key K v (z) is §[27r/3]-nonvanishing. Theorem 3.2(a) 
(with a = 7r/3) implies that Z(G, A, u;x) is §[7r/6]-nonvanishing, and 
Corollary 3.3(a) implies that Z{G, A, u; y 1 ^ 2 !) is §[7r/3]-nonvanishing. 



Example 4.5. Assume that all edge weights are nonnegative reals, 
and that two functions /, g : V — > N are given such that f(v) < g{v) < 
f(v) + 3 for each v G V. Fix the vertex activities as in equation (4.1). 
If every vertex of G has degree at most A then there is a small angle 
e > such that each key K v (z) is §[tt/2 + e]-nonvanishing. To see this, 
the keys with at most three terms pose no problems (by Examples 4.2 
and 4.4). A key with four terms has the form 



J J V/ + V V/ + 2; u + 3 

and the inequality 

d \ / d \ ( d\ ( d 

> 



f + lj\f + 2j \fj \f + 3 

ensures that the only zero of K(z) with nonnegative real part is at the 
origin. Since A is fixed, only finitely many key polynomials need to 
be considered - taking the smallest positive argument of the (nonzero) 
zeros of these to be 7r/2 + e gives the desired angle. 

Theorem 3.2(a) implies that Z(G, A, u;x) is S[e]-nonvanishing, and 
Corollary 3.3(a) implies that Z(G, A, u; y 1 ^ 2 !) is S[2£]-nonvanishing. 

Example 4.6. In Examples 4.1 and 4.2 we concluded that the poly- 
nomial Z(y) = Z(G, A, u; |/ 1//2 1) had only real (and nonpositive) zeros. 
Let Nj = Nj (G, A, u) be the coefficient of yi in this polynomial. It is 
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"folklore" that if Z(y) is S[27r/3]-nonvanishing, then 

NiN k ^ implies that Nj ^ for all i<j<k 

and 

Nf > Nj+iNj^ for all j. 

This property (logarithmic concavity with no internal zeros) is very 
useful for obtaining good approximations to the sequence (Nj) (see 
Eli], for example). 

If all the keys K v (z) are S[57r/6]-nonvanishing then Z(y) is S[27r/3]- 
nonvanishing. However, this hypothesis on the keys is unreasonably 
strong. Consider a key of the form 

with 1 < j < d — 1, corresponding to three consecutive permissible 
degrees. A short calculation shows that this is §[57r/6]-nonvanishing if 
and only if 2j(d — j) < d + 2. This happens only for the pairs (j, d) 
with d < 4 and j = 1 or j = d — 1. 

Nonetheless, I venture the following conjecture. 

Conjecture 4.7. Let G = (V, E) be a finite graph, and let f , g : V — » N 

be any two functions. Fix the vertex activities u as in (4.1). Then 
the sequence of coefficients (Nj) of Z(G, 1, u; y 1 ^ 2 !) is logarithmically 
concave with no internal zeros. 

Example 4.8. Assume that all the edge weights have unit modulus, 
that G is 2/c-regular, and that the key at each vertex is 

(4.2) K(z) = l+^y k + uz 2k . 

If Au > ( 2 ^) 2 then every zero of K(z) has modulus k = u~ 1 ^ 2k . Parts (b) 
and (c) of Theorem 3.2 imply that Z(G, A, u; x) is both kT>- and ac£- 
nonvanishing. Corollary 3.3 implies that every zero of Z(G, A, u; y 1 ^ 2 !) 
has modulus u~ x l k . 

Example 4.9. Assume that all the edge weights have unit modulus, 
and that deg K v (z) = deg(G,v) and every zero of K v (z) has unit mod- 
ulus, for each vertex »gV. Parts (b) and (c) of Theorem 3.2 imply 
that Z(G, A, u; x) is both T>- and £-nonvanishing. Corollary 3.3 im- 
plies that every zero of Z(G, A, u; y 1 ^ 2 !) has unit modulus. 

In particular, these hypotheses evidently hold if A = 1 and the key 
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polynomials are given by K v (z) = 1 + z + z 2 + • • • + z Acg<yG ' v ^ for each 
v £ V. Thus we conclude that every zero of 



E 



If 



FT (deg(G,v)\ 
1 lu€V Vdeg(fi>y 

has unit modulus. 



5. Analogy with statistical mechanics. 

We conclude with an interpretation of Z(G, A, u; y 1 ^ 2 !) inspired by 
analogy with the (canonical ensemble) partition functions in statistical 
mechanics. For simplicity, we restrict attention to a graph G = {V, E) 
that is ci-regular, in which the edge weights A = 1 are all one and the 
activities are the same at every vertex (that is, all the key polynomials 
are equal). The extension to the general case is straightforward. 

The "configuration space" is the set of all spanning subgraphs of G. 
The energy U (H) of a spanning subgraph H C E depends on d + 2 
real parameters J and fi = (/x , fJ>i, fJ>d), as follows: 

d 

(5.1) U(H) =J-#H + Y i iJL j . #Vj(H), 

3=0 

in which Vj(H) is the set of vertices of degree j in H. The quasi- 
physical interpretation of this is that J is the energy of a single edge, 
and fij is the "chemical potential" energy of a vertex of degree j. With 
T > denoting absolute temperature, and (3 = l/k^T where ks is 
Boltzmann's constant, the Boltzmann weight of H is 



(5.2) Z G ((3,J,n) = J2 e 



and the partition function is 

- Will) 
HCE 

This can be interpreted as defining a family of probability measures 
(parameterized by f3, J, and fx) on the set of all spanning subgraphs of 
G: a spanning subgraph H C E is chosen at random with probability 
e~ f3U ^ H ' > /Zc((3, J, /*). A short computation shows that, for H chosen 
according to this distribution, the expected number of edges is 

(5.3) ( #i j) = _^!_iogZ G (/3,J,/x) 
and the expected number of vertices of degree j is 

(5.4) (Wj(H)) = -iAiogZoGM,/*) 
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To continue with the analogy we consider a sequence of graphs G±, 
G*2,... that converges to an infinite, locally finite, limit graph T. (The 
precise definition of convergence is not important for this discussion - 
the prototypical example is that, as n — > oo, the Cartesian product 
C r n of r cycles of length n should converge to the infinite graph Z r 
with edges of Euclidean length one.) We will further assume that the 
"thermodynamic limit" (Helmholtz) free energy 

(5.5) fv(J3,J,^) = ~ lim ) - log Z Gn ((3, J, fj,) 

exists. As in the Lee- Yang theory [SHU], points in the parameter space 
at which the free energy fails to be analytic can be interpreted as phase 
transitions between differing qualitative properties of a random span- 
ning subgraph of T. From the form of (5.5) we see that fr can fail to 
be analytic only at an accumulation point of the union of the zero-sets 
of all the Z Gn (j3, J, /x) (n > 1). From the probabilistic interpretation, 
we are most interested in such accumulation points for which all the 
parameters J, fi) are real. 

The partition functions can be expressed as polynomials in the vari- 
ables 

(5.6) y = e~ /3J and u d = (0 < j < d). 
In fact, a tiny calculation shows that in these variables 

(5.7) Z G (P,J, f i) = Z(G,l,u;y 1 / 2 l) 

with the RHS as defined in (1.4). The point y = 1 corresponds to 
f3J = 0, which is the infinite-temperature limit. If J > then y = 
is the zero-temperature limit, and if J < then y — > +oo is the 
zero-temperature limit. The positive real axis is thus the "physically" 
relevant part of the complex y-plane. If all the chemical potentials fij 
are real then all the activities Uj are positive reals. A zero activity 
Uj = corresponds to an infinite chemical potential fij = +oo, which 
means that a vertex of degree j is forbidden. Notice that the activity 
Uj = e - ^- 7 also depends on temperature except when fij is +oo or 0: 
this is the case precisely when Uj G {0, 1}. 

In this context, Corollary 3.3 has the following immediate conse- 
quence, the proof of which is omitted. 

Proposition 5.1. Let (G n : n > 1) be a sequence of d -regular graphs, 
and let (3 > and J G R and n G M, d+1 be such that the limit (5.5) 
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UjZ J 



exists. Form the key polynomial 

d 

3=0 

with (v,j) as in (5.6). 

(a) If there exists e > such that K(z) is S[ir /2 + e]-nonvanishing then 
/ r is analytic at ((3, J, fx) for all JeR. 

(b) If k > zs swc/i t/iai K(z) is kD -nonvanishing then fr is analytic 
at (j3, J, n) for all 

2 

J > --logK. 

(c) If k > is such that K(z) is kE -nonvanishing and of degree d then 
fr is analytic at ((3, J, fi) for all 

2 

J < --logK. 

Finally, we revisit some of the examples of Section 4, maintaining as 
well the assumptions of Proposition 5.1. 

Example 5.2. With the key polynomial K(z) as in Example 4.3, let 
u = e"^. If ii < +oo (that is, if u > 0) then K(z) is §[tt/2 + e]- 
nonvanishing for all (3 > 0, so that /r is analytic at ((3, J, fx) for all J £ 
R. In this case there is no phase transition at any nonzero temperature. 
On the other hand, if \x = +oo (that is, if u = 0) then both zeros of 

if (2;) have modulus k — ( 2 ) , so that /r is analytic at (/3, J, /*) for 
all 

J>1 f og (t 

In this case there is no phase transition provided that the temperature 
T is sufficiently low compared to the edge energy J. 

Example 5.3. With the key polynomial K(z) as in Example 4.5, 
K (z) is §[7r/2 + e\ -nonvanishing for all (3 > 0, so that / r is analytic at 
(f3, J, n) for all J G K. Thus, there is never a phase transition in this 
model. 

Example 5.4. With the key polynomial K (z) as in Example 4.8, let 
u = and d = 2k. If ( 2 fe fe ) 2 < Au then all the zeros of K(z) have 
modulus k = u~ l l 2k and if (z) has degree d. The only point on the 
positive y-axis at which / r could fail to be analytic is at y = u~ x l k . In 
terms of the "physical" parameters, this says that if 

(5.8) -/3//>21og^ - log 4 
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then a phase transition can occur only at J = —fi/k. The inequal- 
ity (5.8) requires that /i < (so that vertices of degree 2k in H are 
energetically favoured) and that (3 is sufficiently large (so that the tem- 
perature is sufficiently low). If this is the case then a phase transition 
can occur only when the edge energy J and chemical potential /i are 
tuned to satisfy J = —fi/k. 

Example 5.5. With the key polynomial K(z) as in Example 4.9, 
all the zeros of K(z) have modulus one and K(z) has degree d. The 
only point on the positive y-axis at which / r could fail to be analytic 
is at y — 1. In terms of the "physical" parameters, this says that a 
phase transition can occur only at (3 J = - that is, only in the infinite 
temperature limit. 

As these examples illustrate, Proposition 5.1 sees very little about 
the limit graph T - in fact, only the degree of T is relevant. (On the 
other hand, the existence of the limit fr does depend on the structure 
of T.) Thus, for example, Proposition 5.1 can not tell the difference 
between the 3d cubical lattice and the 2d triangular grid - both graphs 
are regular of degree six. Of course, in truth one expects that for any 
given model, the free energies of these two graphs will have different 
phase diagrams. Accounting for more detailed structural properties of 
T remains an interesting open problem. 
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